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Abstract 

The Skyrme-Faddeev model is a (3 + l)-dimensional model which has knotted, 
string-like, soliton solutions. In this paper we investigate a Skyrme-Faddeev model 
with an SO (3) symmetry breaking potential. We then rescale this model and take the 
mass to infinity. This infinite mass model is found to have compact knotted solutions. 
In all of the investigated massive models we find similar charged solutions as in the 
usual, m = 0, model. We also find that their energies follow a similar E ~ Q^/^ power 
growth as the m = model. 



1 Introduction. 



The Skyrme-Faddeev model is a (3 + l)-dimension modified 0(3)-sigma model [Ij. This 
theory has finite-energy string-like solutions and has a potential application as a low energy 
eflFective theory of QCD [2]. The theory usually has two components, a sigma term which is 
a quadratic of derivatives and a Skyrme term which is quartic in derivatives. These are the 
minimum terms required to provide a scale for Derricks theorem, [3j, in (3 + l)-dimensional 
fiat space. The theory is described by the Lagrangian 

L = I (^a^c/> ■ d^cl> - ^(a^c/) X a.c/)) ■ {d^cl> X a'^c/,) - 2mWm^ d'x, (1.1) 

where is a three component unit vector, = (0i, 02,03)- At fixed time is the map 
(p : M.^ ^ S'^. A condition for finite energy solutions of ( |l.lD is that the field must tend to 
a constant value at spatial infinity, which we select to be 0(t, r = oc) = (0, 0, 1) = 63 for all 
time, t. This allows a one-point compactification of the domain, U {0} ^ S^. So static, 
finite-energy, solutions are a map 

0:5^^ S\ (1.2) 
This map, ( 1.2), belongs to an equivalence class characterised by the homotopy group 7r3(5'^) = 



Z. This shows that there is an integer topological invariant associated with 0, known as the 
Hopf invariant. In this case we will refer to the Hopf invariant as the topological charge Q. 
The topological charge can be found by firstly defining an area form u on S^. Then we can 
define g = (j^^uj^ which is the pull-back of uj by 0. Due to the second cohomology group of 
the three sphere being trivial, H'^{S^) = 0, all closed forms on are exact forms. Therefore, 
we can now re-define the exact form g by the 1-form a SlS g = da. Therefore we can express 
the topological charge, as 

Q = ^ [ ^Aa. (1.3) 



4:71 



This non-local definition of the topological charge, (1.3), is not particularly useful in this 



context. Instead the topological charge of (1.2) can be found as the linking number of loops 



in the domain. These loops are formed as preimages of two distinct points on the target 
space. For example if we define the two curves Cp and Cq as the preimages of the points p 
and q. If we then choose a smooth surface D, with a boundary (7^, the linking is defined as 

lmk(C„C,)= J2 ±1' 

DnCq 

where the ± refers to the relative orientations of D and Cq. This definition of the topological 



invariant can be shown to be the same as (1.3), 



There is a well known lower bound on the energy, [SJ [6l [71 18] , which is based on an involved 
argument using Sobolev-type inequalities, 

/ 3 X 3/8 

E>c Q^l^ where c = ( — j • (1.4) 
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The fractional power of (1.4) is proven to be optimal, but the value of c might not be. Ward, 
[6], was motivated by his study of the Skyrme-Faddeev theory on a unit three-sphere to 
propose that c = 1 might be a more optimal value. 

The Hopfion location is commonly defined as the curve C = 0-^(0, 0, -1) = 0"^e_3, which 
is the antipodal value to the boundary vacuum value. 

For m = there have been many extensive and detailed investigations into the static mini- 
mum energy solutions of (1.1), [H ISJ [IHl [HI [III [13]. For charges one to seven it is believed 
their respective global minimum energy solutions have all been identified. It is known [14] 
for topological charges one to three that the minimum energy solutions have a planar loop 
location curve. Topological charges five and six have the minimum energy solutions of two 
linked Hopfions [14]. For topological charge seven the minimum energy location curve is 
a trefoil knot [T4]. In [14j Sutcliffe devised a new knotted rational map ansatz as initial 
conditions. These initial conditions were then energetically minimised to give new minimum 
energy solution candidates for a large class of topological charges. We shall describe and 
make use of this technique later. 



2 m > 0. 



The actual model of interest here is a modification of the usual Skyrme-Faddeev model. It 



is modified by an additional potential term, so m > in (1.1) and 

V((/)) = (l-03). 



(2.1) 



This choice of potential, ( |2.1 ), is not general but is one of the simplest choices. A potential 
term, which breaks global SO (3) symmetry, is also generated on derivation of the theory 
(1.1) as an infrared limit of (3 + l)-dimensional SU{2) Yang-Mills theory [15]. If we restrict 



ourself to the plane this model reduces to the old Baby Skyrmion model. Also this potential, 
(2.1), meets the finite energy criteria for the chosen boundary condition; where the single 
vacuum of V{(/)) is also the chosen boundary value. The potential, (2.1), increases the 



energy density along the location of the Hopfion. This is because the location curve, C = 
0~^(O,O, —1), corresponds to the maximum of the potential (2.1). Therefore, we expect the 
Hopfion location curve to become smaller in length for increasing m. We also expect the 
Hopfion string to become finer with increasing m; as it is analogous with the Baby Skyrmion 
model. This is best understood by an asymptotic analysis where we approximate the field 
for large r as 

where a G (0, 1). For large r we know that, due to finite energy criteria, l^^e^l < 1 therefore 
the energy density associated with ( |l.lD becomes 



Where £„ is a solution of the partial differential equation 

(A - 2m^)ea = 0. 



(2.2) 
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Separating into radial and angular components as = r ^ Ra{r)Qa{d^'4^)- Where Ra{r) 
is a purely radial function, is a spherical harmonic {Vg'^Qa = — A(A + l)©a) and (r, 9^ ip) 



are the usual spherical polar coordinates. Then solving (2.2) we gain an asymptotic approx- 
imation for Ra{r >> 1), 

Where Ca is a constant of integration. This shows that Hopfions become increasingly expo- 
nentially located as m increases. Hence the Hopfions in this massive theory will have Yukawa 
type asymptotic tails. 

3 Initial conditions. 

It is already well known that charged Hopfions can be knotted objects, [H [9l [101 [HI [121 [El 
[T4] . One of the most eflFective ways to create non-trivial knotted initial conditions, to be 
minimised, is to use the rational map ansatz technique described in [T4]. Here the author 
used a degree one spherically equi variant map to compactifiy ^li,Zo ^ 

(Zi,Zo) = sm/,cos/ + — sm/ , (3.1) 

\ r T J 

where /(O) tt, /(oc) and 

-^li^zo = ^o) e C2 I + |Zop = 1}. (3.2) 

Using these complex coordinates an (a, 6)-torus knot can be described as the intersection of 
a complex curve Zq) with a unit three-sphere [1^. Hence we can formulate the rational 
map ansatz 

The inverse stereographic projection of the curve g = produces a 03 = —1 closed curve 



in . The asymptotic value of /(Zi, Zq) in the rational map ansatz (3.3) is used to fix the 
boundary conditions of 0. Therefore, we need /(r)|^^oo = so the inverse stereographic 
projection gives = (0, 0, 1) at the spacial boundary of M^. 



We can now formulate an axially symmetric Hopfion initial condition as 

which has charge Q = np. These types of unlinked Hopfion are labelled An^p. Also a Hopfion 
with a position curve of an (a, 6)-torus knot, can be formed by the rational map ansatz 



W=^, (3.4) 



w = ^° , (3.5) 
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where a G {x > 0|x G Z} and ;5 G {x ^ 0|x G Z}. This gives a closed curve that wraps 
a and b times about the two circumferences of a torus, and has topological charge 

Q = ab + /3a fT4]. These tours knot configurations are labelled A^a,6- A torus knot with a 
reducible denominator can produce linked Hopfions. For example a rational map 

w = ^ = ^ + ^ rs 6) 

Zl-Zl 2(Zi-Zo)^2(Zi + Zo)' ^ ' 



forms two charged n Hopfions linked once. We label this type of configuration as L 



1,1 

n,n' 



This approach not only produces a non-trivial knotted location curve with an analytically 
known topological charge, it also gives a smooth field with the correct boundary conditions. 



To find static Hopf solitons we set = in the Skyrme term of (1.1). This gives a non- 



relativistic theory which has the same static equations of motions as the those derived from 



(1.1). This greatly simplifies the corresponding equations of motion by removing a numeri- 
cally cumbersome matrix inversion. It also still facilitates time evolution by the second order 
dynamics derived from the sigma term. The non-relativist ic equation of motion can then 
be numerically evolved on a discrete lattice using a fourth-order derivative approximation. 
We also need an additional Lagrangian multiplier. A, to constrain to take its value on S'^. 
If we periodically remove kinetic energy the potential energy will also become minimised; 
this will yield minimum energy static solutions. This minimisation technique produces static 
solutions and uses much less CPU time when compared to other similar minimisation al- 
gorithms. A numerical grid of 250 x 250 x 250 points, with Ax = 0.08, was found to be 
large enough to contain the exponentially located Hopfion. On this lattice the Hopfion can 
smoothly attain the vacuum value at the boundary without a noticeable expense of energy. 
This lattice is also fine enough not to lose the fundamental topology. The definition of the 
position of a Hopf soliton is sensible, but not useful for display purposes. Therefore, all the 
images of Hopf solitons in this paper are plots of an isosurface of the preimage of the curve 
03 = —0.85 in the domain. This gives a surface that is a fine tube in the physical space 
and produces much clearer images. To show the linking number we also need to plot the 
preimage of a second loop, but there is no unique loop to choose. In all the plots shown 
we generate general loops on the target space by choosing a point, g = (^1 — jj?^ 0, /i), on 
S'^. We then find the distance on the surface of S'^ between g and 0. This distance is then 
normalised by the distance between g and the south pole of S'^, e_3, as 

COS~-^(g • 0) 



dist(g,0) = 



cos-i(g • e_3)' 



This gives many loops of constant radius on S'^ . The preimages, of these loops, are tubes of 
varying thickness in the domain. Also, an isosurface of unitary value is known to intersect 
with the position curve of a Hopf soliton. Throughout this paper we choose ji = —0.9 which 
is an arbitrary choice for aesthetics. 
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4 Q ^ 4 trivial knots. 



A Hopfion initial condition, where the position curve is contained completely on a plane, 
can be formed by the rational map ansatz (3.4). So if we set n = p = 1 in (3.4) this gives 



the initial condition of a topological charge one Hopf soliton which is located along a planar 
loop. Using this initial condition, and the above minimisation procedure, we can find the 
minimum energy configurations for m = (0, 1, 2, 4, 5), as shown in Figure [l] 




m=0 m = l m=2 m=4 nn=5 

Figure 1: Minimum energy topological charge one Hopfions, all on the same scale. 



CO 




Figure 2: 03 for masses m = 0, 1, 2, 4, 5. 
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Mass, m 


E 





1.236 


1 


1.421 


2 


1.668 


4 


2.017 


5 


2.170 



Table 1: Minimum energies for m = (0, 1, 2, 4, 5) topological charge one Hopfions. 



Figure [T] shows as expected that the larger m is the smaller the Hopfion is. This is due to 
the potential term creating an energy penalty where = (0, 0, —1), which corresponds to the 
maximum value for the potential. Thus the minimum energy Hopfion location loop becomes 
smaller with growing m. The energy for m = is 1.236, this is within 2% of the previously 
accepted result [14j. This similarity is a nice validation of our numerical procedure and the 
small diflFerence can be attributed to the diflFerent choice of lattice spacing. As shown in 
Figure [2] for larger values of mass the field 03 is increasingly symmetric about the Hopfion 
location. Due to this we have decided to perform the remaining analysis with the relatively 
large m = 5. This choice of mass is arbitrary, we could have chosen m > 5, but this choice 
gives sharply located Hopfions which are still large enough so the topology is not lost by the 
numerical lattice. We find the minimum energy Q = 1 Hopfion to be £"1 = 2.17. 
Using (3.4) with (n^p) = (2,1) this again gives a planar loop Hopfion location, but with 
topological charge Q = 2. Minimising this configuration we find the minimum energy 
E2 = 3.45, with a planar Hopfion location curve. We can construct a Q = 3 Hopfion 
by setting (n = 3,]9 = 1) in (3.4) or (a = /3 = 6 = 1, a = 2) in (3.5). The latter configuration 
produces a Hopfion which is located along an unknotted twisted loop. Minimisation of these 
two configurations produces similar configurations with similar energy, = 4.74, located 
on a twisted unknotted loop. 



A topological charge (5 = 4 Hopfion initial condition can be made using (3.5) with either 
(a = a = 2, /3 = 6 = 1) or (a = a = 4, 6 = 1, /3 = 0). This gives Hopfion location curves 
K2^i and ^4^1 respectively. Minimising both of these configurations give a ^4^1 (twisted ^4,1) 
Hopfion location curve, with energy E4 = 6.051. We can also generate Q = 4: axial initial 
conditions using (3.4) with {n = p = 2) and (n = 4,]9 = 1). This gives ^42,2 and A^^i planar 
curves respectively. Minimisation of the ^2,2 configuration remains as an ^2,2, with an energy 
^ 0.6% larger than the ^4,1. This is within numerical accuracy of our minimisation scheme. 
Therefore we are not able to define which of these two configuration is the lower energy. The 
^44^1 configuration also minimises to a planar curve described by A^^i. This seems to show 
that twisting the loop slightly reduces the energy. The Q = 4: planar loops are most likely 
long lived saddle point solutions, preserved by symmetry. 



7 



oOO 



Q=l 0=2 0=3 0=4 

Figure 3: Minimum energies for m = 5, Hopfions. 

5 Q ^ 5 knotted/ linked Hopfions. 

For a select few simulations Figure [4] and [5] show a plot of minimum energy Hopfion locations. 
Also Figure [4] and [5] show the linking of the initial rational map ansatz and the linking of the 
corresponding minimum energy Hopfion curve. For topological charges (5^4, both in the 
massive case and in the normal massless case, the Hopfion location curves are all found to 
be unknot solutions. For topological charges Q ^ 5 we find the minimum energy Hopfions 
have either a linked or knotted location curves, as shown in Figure [5| The minimum energy 
solutions presented have very similar qualitative features with the massless model [14j. The 
minimum energy Hopfion location curves for each charge sector have comparable linking form 
hence, due to the computational intensity, we have restricted our analysis to the presented 
charges. 
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o 



7^5,2 i^3,2 7^^4,3 ^ i^2,3 7^7,1 ^ ^7,1 





€) O 



5L^;; ^ L'^l bKs,2 ^ lJ;^ 5i^:4,i ^ ^5,1 



Figure 4: Position curves for m = 5 Hopfions of charge 5 < Q < 7. Showing their charge 
and hnking form. 
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^ ^ ^4,3 16^4,3 ^ i^4,3 




Figure 5: Position curves for m = 5 Hopfions of charge Q = (8, 12, 15, 16). Showing their 
charge and hnking form. 
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E 




Q = 4 






K2,i 44.1 


6.05 


2.14 




6.05 


2.14 


A4,l ^4,1 


6.07 


2.15 


A.2,2 ~^ ^2,2 


6.09 


2.15 


7-1,1 ,-1,1 


7.17 


2.53 


Q = 5 






7-1,1 .1,1 
-^2,1 ^ -^1,2 


6.23 


1.86 


-f'^3,2 


7.17 


2.14 


A4I ^ ^5,1 


7.43 


2.22 


Q = 6 






-f^3,2 — > ^3,2 


8.01 


2.09 


^2,2 L^'o 


8.19 


2.14 


7-1,1 .1,1 
-^3,1 ^ -^3,1 


8.41 


2.19 


^5,1 ^ i^5,l 


8.60 


2.24 


^4,2 ^ i^Ji 


9.07 


2.37 


Q = 7 






i^5,2 ^ i^3,2 


9.19 


2.14 


-^3,2 ^ -^3,2 


9.20 


2.14 


^4,3 K2,3 


9.64 


2.24 


^7,1 ^ ^7,1 


10.18 


2.37 


g = 8 






-^2^2 ~^ ^4,2 


9.86 


2.07 


-f^3,2 — > ^4,2 


9.88 


2.08 


i^3,4 ^ ^4,2 


9.88 


2.08 


^5,2 ^ ^4,2 


9.97 


2.10 


rl,l rl,l 
-^3,3 ^ ^3,3 


10.45 


2.20 


Q = 12 






-f^4,3 1^2^2,2 


13.77 


2.16 


-f^5,3 ^ -^2'2'2 


13.77 


2.16 


7-1,1 ,-1,1 


15.30 


2.37 


Q = 15 








16.17 


2.12 


/S:4,3 i^4,3 


16.56 


2.17 


-f^3,2 ^ -^7^2 


17.11 


2.25 


Q = 16^ 






7-3,3,3,3 7> 

-^1,1,1,1 ^ -'^4,3 


17.07 


2.13 


i^4,3 ^ i^4,3 


17.12 


2.14 


^3,2 K4,3 


17.25 


2.17 



Table 2: m = 5 Hopfions initial and final configurations, with their respective energies. 



11 



One main difference between our results and those of [T4] is that we are presenting po- 
tentiaUy new energetic local minimum or saddle point solutions. As shown in Figure [5] we 
have new linked topological charges Q = 8, 12, 15 excited solutions. 

Also our topological charge eight and six minimum energy solutions seem to have similar 
linking structure as in the m = case. But in the m = 5 case the links seem to be almost on 
top of each other. The minimum energy location curve of the topological charge 15 Hopfion 
is topologically similar to that found in [Hj, but it is qualitatively different. Figure [T] shows 
that the energy, as a function of topological charge, has a similar growth to the m = case 



6 Infinite mass, m oo. 



Rescaling the Lagrangian density (1.1) by ^ ^/mx gives 



C{y/rnx) C2 ^ ^ 

1= = h + vLq. 

'm m 



Where Ca refers to the a^^ order derivative in the Lagrangian density. We can then define 



lim — — = £4 + £o- 



m 



(6.1) 



A model comprising only of a Skyrme term and a potential term has been addressed before 
[TTj . In [TTj this model was derived by setting a constant to zero. This effectively removes the 
sigma term in the Skyrme-Faddeev with a potential model. Soliton solutions of this model 
are commonly referred to as compactons [TTj . This is because they reach their vacuum value 
in finite distance and therefore have no asymptotic tails. Hence they are effectively BPS 
for large separation. In order to numerically find minimum energy solutions of (6.1) it is 
computationally easier to work with a modified model of the form 



do<p • dpcj) 
327r2^ 



+ A 



1 9*0=0- 



(6.2) 



This modified model can be simulated by a trivial extension of the previous Skyrme-Faddeev 



model. The equation of motion and the energy density of (6.2) will converge with the infinite 



mass case, (6.1), in the static limit. This model, (6.2), also allows for time evolution by the 



second order dynamics of the purely kinetic sigma term. Again we use the rational map 



ansatz, (3.5), to give non-trivial knotted initial conditions. For this infinite mass compacton 



model, (6.1), we find the topological charge-specific minimum energy candidates in Table ^ 



The numerical scheme is fundamentally the same as the one used in the finite mass case. 
But now due to this model having a different scale we found Ax = 0.1 to be a suitable lattice 
spacing. The results and initial conditions of this investigation are shown in Table [3} Also, 
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for topological charge 1 ^ Q ^ 5 the Hopfion location curves are shown in Figure (H) We 
found that the minimum energy solutions for 1 ^ Q ^ 5 are similar to the m = 5 and m = 
models. The compact nature of this model can be seen in Figure |8} which show 03 along a 
radius that is in the same plane as the planar Q = 1 Hopfion. This shows how the compact 
Hopfion field attains the vacuum value in finite distance. This shows that in the m = oo 
model two well-separated static Hopfions do not attract or repel each other. This is due 
to there being no overlap of the Hopfion tails. Therefore, the string self-interaction of this 
m = oo model is much less than in the finite mass model. Also the shape of 03, in Figure [8| 
is more symmetric about its minimum than the plots of 03 in the finite mass model. Figure 
[2| This shows that in the m = oo limit the string cross-section is much more symmetric. 
Again, in this infinite mass case, we have found new local minimum energy solutions for 
(5 = 4; both of which have location curves described as a very twisted ^44^1. A point worth 
noting is that for these solutions the boundary of the numerical lattice is very far from the 
Hopfion. 





P 
j-j 




Q = l 








0.86 


0.86 


Q = 2 






A2,i A2,l 


1.37 


0.82 


Q = 3^ 






K2,l is,! 


1.90 


0.83 


^3,1 ^3,1 


2.02 


0.90 


Q = 4 






^2,2 ^2,2 


2.50 


0.88 




2.56 


0.91 


K2,l ^4,1 


3.19 


1.13 


rl,l . 
-^1,1 ^ -^1,1 


3.38 


1.20 (not shown) 


K4,l ^4,1 


4.22 


1.50 


Q = 5 






7-1,1 ^1,1 

-^2,1 ^ ^2,1 


3.00 


0.90 




3.23 


0.97 


^3,2 — ^ -f^3,2 


3.88 


1.16 


Q = 7 






k3^2 ^3,2 


3.97 


0.92 



Table 3: Infinite mass Hopfion initial and final configurations, with their respective energies. 
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5^5,1 ^5,1 





4^2,2 ^2,2 



4^4,1 ^ ^4,1 



4i^2,l ^ 





3K2A ^ ^3,1 




3^3,1 ^ ^3,1 





2^2a ^ ^2,1 



IA1 1 ^ ^11 



Figure 6: Position curves for m = oc Hopfions of charge 1 ^ Q ^ 5. Showing their charge 
and hnking form. 

Table [3] shows that in this m = oc model the initial and final minimised location curves 
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Figure 7: The ratio of the energy to the topological charge one Hopfion, as a function of 
topological charge Q for the three models. 



rarely differ. A good example of this is the Q = 5 K^ 2 trefoil knot which under minimization 
remains as a trefoil knot, but in the m = 5 and m = models minimises to a L^^^^. 
The (5 = 5, L2'i, is also a lower energy solution for this m = oo model. This reluctance 
to deform from one location curve to another is due to the reduced self interaction of this 
compact Hopfion model. 
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1 2 3 4 5 

r 

Figure 8: 03 along a radius of a Q = 1 Hopfion, in the infinite mass case. 

We have also found a topological charge seven trefoil knot, in the m = oo model, which 
is shown in Figure [9| For the range of topological charges we have investigated, in this 
m = oo model, we have found that the charge-specific lowest energy solutions are similar to 
the m = and m = 5 models. But as in the m = 5 model we have discovered more local 
minima or saddle point solutions. Again, as shown in Figure [7| the energy per-unit charge 
scales with a similar E ^ Q^l^ grows as the m = [H] and m = 5 models. 
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Figure 9: The Q = 1 trefoil knot minimum energy solution, for the m = oo model. 

7 Concluding remarks. 

We have explored the Skyrme-Faddeev model with a potential term included. We found 
that including a potential in the model makes the Hopf solitons exponentially decay to their 
vacuum value. Increasing the mass makes the Hopf solitons string cross-section increasingly 
exponentially localized. We found for m = 5 the minimum energy solutions are described 
by similar linking curves as the m = Skyrme-Faddeev model [14j. 

Using a spatial rescaling we were able to formulate an infinite mass model. This infinite mass 
model is known to yield compact Hopfions [17j. For this infinite mass compact model we have 
presented a number of topologically charged solutions. We showed that the minimum en- 
ergy compact solutions have similar location curves as in the usual, m = 0, Skyrme-Faddeev 
model and the m = 5 massive model. In both the m = 5 and m = oo models we found 
new local minimum, or stationary point energetic solutions. The increasing localization of 
the strings in the two massive models reduces the string self interaction. Therefore, it is not 
surprising that these models possess solutions stabilised by symmetry. 

In Figure [7| we showed that the m = 5 and the m = oo models charge-specific energies 
seem to grow with a similar E ^ Q^^^ trend as the m = case. Therefore, approximating 
Hopfions as KirchhoflF rods [18j could also be successful in the massive model. 
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